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Abstract 

We use the recently developed massive field approach to calculate the scalar self-force on a static 
particle in a Schwarzschild spacetime. In this approach the scalar self-force is obtained from the 
difference between the (massless) scalar field, and an auxiliary massive scalar field combined with 
a certain limiting process. By applying this approach to a static particle in Schwarzschild we show 
that the scalar self-force vanishes in this case. This result conforms with a previous analysis by 
Wiseman |9|. 
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I. INTRODUCTION 



The self-force is a force originating from the coupling between the charge of a particle and 
the field that this charge induces. In curved spacetime, one usually considers a fixed curved 
background spacetime and a given worldline of a charged particle. The field induced by 
the point particle diverges at the particle's location, and therefore a regularization method 
is required to calculate the correct (and finite) self-force from this diverging field. For this 
problem, formal expressions for the self-force have been derived for various types of fields. 
Thus, DeWitt and Brehme derived an expression for the electromagnetic self-force 111. More 

□ n 

recently Mino, Sasaki and Tanaka |2J, and independently Quinn and Wald [3J, have derived 
an expression for the gravitational self-force (here the electric charge is replaced by the mass 
of the particle, and the induced electric field is replaced by the linear gravitational pertur- 
bation induced by the particle's mass). Quinn jj| has recently derived the corresponding 
expression for the scalar self-force. 

Although formal analytical expressions for the various types of self-forces are well known, 
explicit analytical calculations of these self-forces are still a challenging task and were carried 
out only in few cases. Here, one of the main difficulties is the non-local nature of the self- 
force in curved spacetime (i.e., the self- force value at each point on the particle's trajectory 
generically depends on entire past history of the particle). For a weakly curved spacetime 
the above mentioned explicit self-forces expressions were derived by DeWitt and DeWitt 
j^j], and by Pfenning an Poisson 6]. The self- force on a static particle was investigated 
analytically by several authors: Smith and Will have obtained a non-vanishing result for 
the electromagnetic self-force on a static particle in Schwarzschild [7]. Later Lohiya ^] 
derived the electromagnetic self-force on a static particle for other types of background 
spacetimes. Recently, Wiseman t 9] has showed that the scalar self-force on a static particle 
in a Schwarzschild spacetime vanishes. 

A practical method to calculate the self-force for generic orbits was devised by Barack 
and Ori (this method was later improved U,^^). In this method, one first calculates 



certain regularization parameters (usually analytically), and then uses these parameters to 
calculate the self-force (usually numerically). This method was implemented numerically in 



certain cases 



L13, 



14j . In particular, using this method (numerically) the scalar self- force on 
a static particle in Kerr-Newman background was calculated by Burko and Liu For 
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other approaches to the self-force problem see [ly, |17 

Very recently, a new general method titled "the massive field approach" for the calcu- 
lation of the scalar self- force was developed Q]. In this paper, we implement this new 
method to calculate analytically the scalar self-force on a particle, which is held static (by 
some external forces) in a Schwarzschild background spacetime (i.e., static with respect to 
Schwarzschild coordinates), and thereby we show that the self-force vanishes in this case. 
This result conforms with the above mentioned analysis by Wiseman (which used a differ- 
ent calculation method). The analysis given here therefore verifies Wiseman's result and 
furthermore demonstrates how the massive field approach can be implemented for this type 
of self-force calculations. 

In Ref. [19] it is shown that the scalar self-force can be obtained from the difference 
between the following retarded scalar fields: The first field is - the massless scalar field 
induced by the particle. This field satisfies the inhomogeneous massless scalar field equation, 
with a charge density p [see Eq. (jH) below]. The second field is a - an auxiliary massive 
scalar field satisfying the inhomogeneous massive field equation, with the same charge density 
p [see Eq. (j!4j) below]. Ref. [l9j provides the following prescription for the calculation of 
the scalar self-force f^ el ^ in terms of these fields 

f; elf {zo) = q Km (lim A0,» + lq[A 2 n^(z ) + Aa,(z )] ) , (1) 



where, 



Acf>(x) = <f>{x) - A (x) . (2) 



Here, A is the mass of the massive field, Zq is the self-force evaluation point on the particle's 
worldline, and x is a point near the worldline defined as follows: at zq one constructs a 
unit spatial vector n M , which is perpendicular to the particle's worldline but is otherwise 
arbitrary (i.e., at zq the following relations are satisfied: n^n^ = 1, n^u^ = 0). In the 
direction of one constructs a geodesic, which extends out an invariant length 5 to the 
point x(2o,n M ,5). In the following section we shall calculate the scalar self-force on a static 
particle in Schwarzschild spacetime by implementing the mathematical prescription given 
by Eq. JT}. 
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II. IMPLEMENTATION OF THE MASSIVE FIELD APPROACH 

Here we use Eq. (jj) to calculate the scalar self-force on a static particle in 
Schwarzschild. We use Schwarzschild coordinates throughout, where x a = (t,r,6,<p), 
g^ u = diag(— f, r 2 , r 2 sin 2 6), / = 1 — — . We denote the particle's worldline with 
z(t), where r denotes the particle's proper time. More explicitly the particle worldline is 
characterized by the spatial coordinates (r ,#o,y2 ) which are constants, where r > 2M. 

By virtue of the spherical symmetry of the problem the two angular components of the 
self-force fg el ^ and vanish. To calculate the time component of the self-force f* el f 
we use Eq. (|T|). For a static particle in Schwarzschild the scalar fields <fi and a are time 
independent, and therefore the term A0 jt in Eq. ((TJ vanishes. Noting that the vectors 
and a M are perpendicular to the particle's worldline we find that f^ el ^ vanishes, by virtue of 
Eq. (dJ. 

We now consider the radial component of the self-force f^ el ^ . Eq. implies that the 
calculation of f^ el ^ follows from the calculation of the field A0 jr near the particle's worldline, 
combined with a limiting process. The first limit that has to be considered is 5 — > 0. Here, 
we take advantage of the arbitrariness in the definition of the unit spatial vector n M , and 
choose n M be the radial unit vector = (0, /q -1 ^ 2 , 0, 0), where /o = /(?"o)- This choice 
defines the spatial geodesic which extends out from zq to x in the direction of to be a 
radial spatial geodesic (with r > r ) denoted here with 7. The limit S — > is now equivalent 
to the limit r — > r$ along this geodesic. From Eq. Q we now find that the radial component 
of the self-force is given by 

fff(z Q ) = q lim ( lim A0, r (r) + ^[A 2 n r (^ ) + Aa r (z )] ) . (3) 

A^oo I r-*r 2 J 

Here the limit r — > r is considered along 7. 

This section is organized as follows: First, in Sec. Ill Al we calculate the field jr along 
the radial geodesic 7, then in Sec. Ill Bl we calculate an expansion for the field 0A,r i n the 
vicinity of the particle's worldline and along 7, and finally in Sec. Ill CI we substitute 4> >r and 
0A,r in Eq. Q and calculate f^- 
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A. Massless scalar field 



Here we calculate jr along the spatial radial geodesic 7. The massless scalar field 
satisfies 

□0 = -4tt p . (4) 

Here □</> = ;A /\ and p(x) is the scalar charge density. For a point particle this charge 
density is given by 

poo -i 

p(x) = q 6 4 [x - z(r)]dr , (5) 



-0 

where g denotes the determinant of the background metric. For a static worldline Eq. (JSJ 
gives 

qV7o$(r - r o )5(0 - 6» )<5(^ - <p ) 

P = — — a • 6 

r l sin U 

Next, we solve of Eq. (jlj) for the charge density given by Eq. (jHJ), for this purpose we 
decompose p and into spherical harmonics: 

00 1 

Z=0 m=-Z 
00 / 

z=o m=-; 

The coefficients in the charge density decomposition are given by 

Am <H r ~ r o)v / 7o^ m 



P = 1 r ^ WJU Y lm *(6 , <po) . (9) 

The above mode decompositions together with Eq. (jlj) give the following set of decoupled 
ordinary differential equations (for a static source) 

{r 2 f^) r - 1(1 + l)0' m = -4nr 2 p lm . (10) 

These equations have the following analytical solutions 

-Y^Oa^PiMQ^). (11) 



dm \/-lm* 1 



M 

Here z = — 1, 20 = - 2 ( r o); - 2 > an d z < denote the larger and smaller terms from the 
pair {z , z}, respectively; and Pi, Qi denote Legendre functions of the first and second kind, 
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respectively. We now substitute Eq. (fTTj) into Eq. (JHJ) and sum over the multipole number 
m which gives 

^ = E( 2/ + l)Pi(z<)Qi(z>)Pi(™sa) . (12) 



M 

1=0 



Here, 



cos a = cos #0 cos 9 + sin 9 sin 9 cos(y? — y^o) . 



Recall that here we are interested only in the solution along the radial geodesic 7. Along 7 
we have cos a = 1, which considerably simplifies the summation in Eq. ()12j) . and we find 
that along this geodesic 

(p = q- 



r — r 

Finally, we differentiate this expression with respect to r and obtain 



-"V^kf- (13) 



B. Massive scalar field 



We now calculate an approximate expression for 0A,r along the radial geodesic 7. Note 
that for the implementation of the prescription summarized by Eq. (jSJ), it is sufficient to 
have an approximate expression for <fix,r in the vicinity of the particle's worldline, and as 
A — > 00. We therefore expand 0A jr along 7 in powers of (r — ro) and in powers of A" 1 ; and 
keep only terms that do not vanish as r — > r , and as A — > 00. 

The massive scalar field 0a satisfies 

(□-A 2 )0A = -4vrp. (14) 

Here, the charge density p is given by Eq. © - the same charge density as in the massless 
field equation. Decomposing 0a, and p into spherical harmonics gives 

00 1 



E E -y lm (0, V )^(r), (15) 

z=o m=-l 
co I 1 

P = E E -y lm (0^)P lm {r), (16) 



Z=0 m=-l 

where p lm = rp lm . From the massive field equation ()14j) together with the spherical har- 
monics decompositions ()15I16|) . we obtain the following infinite set of decoupled ordinary 



differential equations for the spherical harmonics coefficients <p l A : 



frr* - (V\r) + A 2 fW = ~^p lm f . (17) 



Here we introduced the tortoise coordinate r* [see e.g., j^] defined by ^ = j, we also de- 
fined V l (r) = /(— + -^—5-^) and /' = j-. We now express the solution of the inhomogeneous 
equation ()17|) in terms of two linearly independent solutions of the corresponding homoge- 
neous equation. Denoting these homogeneous solutions with x± : such that x+ vanishes as 
r* — > 00, and x~ is regular as r* — > —00, we find that 



A ' _ r W\xl,X l +}r 



(18) 



Here = X+ *X- ~~ X- *X+ denotes the Wronskian, and r > and r< denote the 

larger and smaller terms from the pair {r ,r}, respectively; the subscript r indicates that 
the term in the square brackets is evaluated at r = tq. 

Next, we calculate the sum over the multipole number m in Eq. (|15|) . Introducing 
<f) l A = Ylm=-i l^f^Y 1 ™ 5 an d using Eq. (fTH|) we find that along the spatial radial geodesic 7 
the multipoles (f) l A are given by 



, = -2gL v ^oX-(^o)x+( r ) 
A "" rr W[xi,X l + U 



(19) 



where L = I + |. 



Next, we calculate an asymptotic expansion for 0A r as A — >• 00. This requires summation 
of the terms <j) l A r , followed by an asymptotic expansion in powers of A -1 . First we consider 
the summation over the multipole number / 



00 

'a.,- =TAa,- (20) 



1=0 

Recall that for our purposes it is sufficient to calculate 4>A,r along the radial geodesic 7, and 
up to 0[(r — r )°] (inclusive). As discussed in Ref. j3| the field 0a,t diverges on the particle's 
worldline. However, the individual multipole terms 4> l \ r {ro), are finite on this worldline 1 , 
only their sum over the multipole number / diverges there. To calculate this sum (off the 
worldline) it is useful to split this sum into two parts: a sum which diverges at the limit 



1 For an approximated expression for these multipoles see Ea. l25|) below. 
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r — > ro, and a sum which is finite at this limit. This splitting simplifies the calculations, and 
allows us to use different calculation methods for the two sums. We therefore express (p\ )r 
as 

oo oo 

0A,r = £>' + -h l ). (21) 

1=0 1=0 

This splitting is considered on the radial geodesic 7, and the functions hi(r) will be defined 
below, such that the second sum in Eq. (}2*T]) remains finite as r — > r$ [see Eq. (13 1|) below]. 
This requirement implies that the values of hi(r) in the vicinity of the worldline must reflect 
the leading asymptotic expansion of 4> l \ r {r) as L —>■ 00. We comment here that a similar 
method for calculating a mode-sum was devised in Ref. [11] (for a somewhat different 
purpose). 

We shall now derive the required expressions for the functions h (r). This requires an 
analysis of the asymptotic behavior of <p l A r (r) as L — * 00 in the vicinity of the worldline. For 
this purpose we use the WKB approximation, which enable us to calculate an asymptotic 
expansion of the solutions of Eq. ()17|) in inverse powers of L (or in inverse powers of A). 
Here there is a problem though, since the WKB approximation is invalid in the vicinity 
of the black-hole horizon. This can easily be inferred from Eq. (fTTjl which has a "turning 
point" at the event horizon [i.e., the term in the brackets in Eq. ([17]) vanishes at r = 2M]. 
Therefore, we can not impose a boundary condition at r = 2M on the (WKB approximated) 
solution We shall therefore consider a general approximated solution for x_ m the region 
where the WKB approximation is valid, without imposing any boundary condition. We 
shall eventually deal with this boundary condition problem in appendix A, by considering a 
different approximation method. For the accuracy required by the calculations below it will 
be sufficient to keep the first three leading terms in the WKB approximation, which reads 

x l + « e ~ 5o+5l ~ 52 , (22) 
X l _ ~ Cie + 5 o+Si+s 2 + C2e -Sb+S!-s 2 _ (23) 

Here c\ and C2 are independent of r*, and the functions Si, (i = 0, 1, 2) are given by [see e.g. 



S 



21]] 



S, 



o 



[' f-'VUdr', (24) 

J2M 



Si = ~lnl7, 

r i f /w _ s(/^) 2 . 

hu! \ ~ 321/6/2 1 " • 



where U(r') = V l (r') + A 2 /(r'), and U' = jy. Examining the asymptotic properties of the 
functions d r Si, we find that as / — > oo (for a fixed A) the functions d r Si are 0(-L 1_ *); and as 
A — > oo (for a fixed L) the functions <9 r Sj are 0(A 1- *). We now substitute equations (|22|2Hj) 
into Eq. p9j) . and differentiate with respect to r, this gives 



Le 5 . 1 



forr [S 0>r + S 2 



(S,r - -) , (25) 



,r\r 



r 



where, 

S = -[S (r) - S (r )] + [S^r) - S 1 (r Q )} - [S 2 (r) - S 2 (r )] . (26) 

This WKB approximation is accurate up to 0(/ _1 ) (and up to 0(A -1 )). At the limit where 
/ — > oo (or A — > oo) the contributions to Eq. (|2*3jl from the term which contains the 
coefficients c\ and c 2 [originating in Eq. (|23|) ] vanishes faster then any negative power of L 
(or A), and therefore this contribution was neglected here (see Appendix A for details). 

To find the asymptotic expansion for </>^ r ,as / — > oo, we first calculate the leading 
asymptotic expansion of S. Expanding S gives 

S = -aL + 0{L°) . (27) 

Employing equations (|24l2b127|) we find that a is approximated by 

Eq. (|27|) implies that as a function of L the term e s in Eq. (J25j) has an essential singularity 
at infinity. To deal with this singularity we simply multiply Eq. (J25j) by e~ aL e aL , and obtain 



i w e~ aL 



Jorr [So tr + S 2) 



rjr 



r 



(29) 



Note that no further approximation was made here, since we merely multiplied Eq. (|2~5j) by 
unity. However, in this form the term in the square brackets does not contain an essential 
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singularity as L — > oo, and therefore this term can be expanded in powers of L without 
difficulties. Note that higher orders in the expansion of S do not give rise to a similar 
essential singularity, and therefore do not require special attention. Moreover, higher orders 
in the expansion of a in Eq. (J28j) will have a vanishing contribution at the limit r — > r (see 
below), and are not required here. An expansion of the term in the square brackets in Eq. 
(JUJ) reads 

q Le s+aL (5 , _1, 

V7o rr O [So,r + S 2, r ]r ^ V 

Here the coefficients A(r),B(r) and C(r) are independent of L. We now define the functions 
h l (r) by multiplying the asymptotic expansion in Eq. IjflOj) by e~ aL , which gives 



A(r)L + B(r) + 



C(r) 



+ 0(L- 



(30) 



h\r) = e~ a{r)L 



A(r)L + B(r) + 



C(r) 



(31) 



Having this definition, the functions h l (r) coincide with the leading asymptotic expansion of 
^Ar( r ) i n ^ ne vicinity of the worldline, as required. Moreover, the difference (<f) l Ar — h l ), when 
evaluated on the worldline is 0(L~ 2 ). Therefore, the second sum in Eq. (|2*Tj) converges on 
the particle's worldline, as required. We comment that at the limit r —* r the coefficients 
A(r), B(r), C(r) coincide with the first three mode-sum regularization parameters introduced 
by Barack and Ori in Ref. nj (for the particular problem which is considered here). 
Next, we calculate the first sum in Eq. (|2T|). Employing Eq. (|3"T|) we obtain 



-aL 



AL + B+ C y 

Li 



(32) 



1=0 1=0 

Note that for r > ro this sum converge due to the factor e~ aL . Summing separately over 
the various terms (without the coefficients) in the square brackets in Eq. (J32|) gives 

aL _ cosh(«/2) 



1=0 
oo 

1=0 

oo 

E^- 1 



4sinh 2 (a/2) 
2 



(33) 



sinh(o;/2) ' 
- aL = 2arctanh(e~ Q/2 ) 



1=0 



Eq. (j2Hj) implies that as r — > r the first sum in Eq. is 0[(r — r ) 2 

sum is 0[{r ! _ 1 1 ■ 1 ' ' ' ' : - 1 : ' 1 1 1 1 ' 1-1 



the second 

/ni "J, and the third sum is 0{log[(r — r ) _i ]}. Our calculation has to be 



10 



accurate up to 0[(r — r )°], and therefore the coefficients A(r), B{r), C{r) in Eq. (J3*2|) . and 
the function a(r) have to be accurate up to [0{r — r ) 2 ]. We now substitute Eq. (p?3"j) in Eq. 
(132)) . and calculate the coefficients A, B and C to the required order. For this calculation 
we use Eqs. (I24l2(il28l3()|) . We find that the first sum in Eq. fl2]} is given by Q 

£>' = -r^W + 77r A2 + °(' - r °) ■ ( 34 ) 

^ (r - r ) 2 2^ 

Next, we consider the second sum in Eq. (|21|) in the vicinity of the worldline: 

oo oo 

5>a>) - ti(r)] = $>i, r (r ) - h l (r )] + 0{r - r„) . (35) 

Recall that we only need the asymptotic form (as A — > oo) of this expression. Introducing 
the notations y l = l/A, y l = L/A, and </>A >r (L, r) = 0^ r (r), h(L,r) = h l (r). We find that as 
A — > oo the sum in Eq. ()35j) can be approximated with a Riemann integral 



^ /"OO 

A^A,r(y r A,r )-/i(|/ r A,ro)](j/ ,+1 -j/ , )«A / [<j>A,Mv,r ) - ti(Ay,r )]dy . (36) 

We now substitute the WKB approximation given by Eq. ({25} into the integral in Eq. (|3fij) . 
Note that contributions to (j)\ jr (y A, Tq) from the functions with « > 3 are convergent 
upon summation over the multipole number I. However, these sums vanish as A — > oo, and 
therefore it is sufficient to keep only the first three leading terms in the WKB approximation. 
We now substitute Eq. ()25j) into Eq. (|3l)j) and expand the result in inverse powers of A 
[technically it is useful to substitute U(r) = U + AU(r) into Eq. (J2HJ), where U = U(r ), 
and formally expand this expression with respect to Uq, and keep only terms that will 
eventually give a non-vanishing contribution as A — > oo] we find that in the vicinity of the 
particle's worldline the second sum in Eq. ()21|) is given by 

qM 

\<PA,r -">') = ; 

1=0 

Substituting equations (|34I37|) in Eq. (|2*T|) gives 

(l\/To , q 



22 



£(0A, r ~ h l ) = + 0(A' 1 ) + 0(r - r ) . (37) 



>A,T 



(r — tq) 2 2 



+ ±(Aa r + AV) + 0{r - r ) + O^ 1 ) . (38) 



Here n r = f is the radial component of the unit spatial vector n M , and the 
radial component of the four-acceleration of the static particle in Schwarzschild. 
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C. The radial self-force 



We now calculate the radial self- force by implementing the prescription given by Eq. Q. 
First we calculate the field A0 r (r) in the vicinity of the particle's worldline. From equations 
(USEHD we obtain 

A0, r = -|(Aa P + AV) + 0(r - r ) + O^ 1 ) . (39) 

Note that the field A0 )T .(r) remains finite as r — > ro (while the fields ir and 0A,r diverge at 
this limit) - this is a general property of A0 r (r) (see Ref. |l9j). The cancellation of the 
divergent terms allows us to take the limit r — > r of A0 r (r). Following Eq. (J3j) we add the 
term |[Aa r (zo) + A 2 n r (^o)] to this limit. This added term exactly annihilates the 0(A) and 
0(A 2 ) terms in our expression - this is a general property as well (see Ref. [3]). We now 
complete this calculation by taking the limit A — > oo in Eq. (J3J), which gives 

/ r sd/ (z ) = g lim ( lim A0, r (r) + \q[k 2 n r {z Q ) + Aa r (^ )] 1=0. (40) 

Combining this result with the results for the other components of the self-force, we conclude 
that the scalar self-force on a static particle in Schwarzschild vanishes. 
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APPENDIX A: GREEN-LIOUVILLE APPROXIMATION 

Here we shall derive Eq. (125)1 . In particular we shall provide justification for neglecting 
the term which depends on c\ and c% in the expression for <fi l Ar . As discussed in subsection 



III B| the WKB approximation is invalid in the vicinity of the event horizon. This situation 
prevents us from patching the WKB approximation for the homogeneous solution x l - to a 
regular boundary condition at the event horizon. In the region where the WKB approxima- 
tion is valid, the solution x l - is given by Eq. (j23j) . By substituting equations (j22|23j) into Eq. 
(|19j) . and differentiating with respect to r we find that 



V/oTo P0,r + ^2,rjr r Ci 



-(SV--)(l + ^e- 2S °- 2S2 ). (Al) 
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We now focus on the contribution from the last term in the last brackets, and explain why 
;his term can be neglected here. For this we follow an analysis by Rowan and Stephenson 



231 ] . who showed that by using the Green-Liouville approximation method, one can obtain 
an approximate solution to Eq. (|17j) . which is valid in the entire r > 2M region. First, we 
express the homogeneous wave equation for x± as 

[r 2 fx l ,rir-W + l)+A 2 r 2 ] X l = 0. (A2) 

Next we introduce x = (r/M)f, N = MA; and make the following transformations from x 
to £: 

-/2 ^9/2 + x x , a 2 



x x(x + 2) 

Here, f = (f ), /3 2 = N 2 /k 2 , a 2 = [1(1 + 1)}/ k 2 , and k 2 = N 2 + l{l + l). Note that < a < 1, 
< (3 < 1. The new variable £ is defined by 



i= [ X \e\dx. 
Jo 



Next, we define ip l such that 



i, 1 = ^'x(x + 2) X l . (A3) 

Eq. ()A2|) now reads 

4 5 -[fc 2 -(2O- 2 + ^i(O]^ = 0. (A4) 

Here ^i(^) is a slowly varying function of £, which is bounded everywhere. This function is 
Off) near the event horizon, and it is 0(£~ 2 ) as £ — > oo. The full expression for gi(£) (see 
231 ]) is not required here. 

For large values of k (which corresponds to large values of N and/or large values of /) 
Eq. (|A4|) can be solved using perturbation analysis. The leading order in this approxima- 
tion is obtained by completely neglecting the small contribution from the function #i(£) 
in this equation. The remaining equation can be solved exactly, giving two linearly inde- 
pendent solutions, from which the corresponding x± can be calculated. These approximate 
homogeneous solutions are given by 

X\ « ViWx{x + 2)}~ 1 / 2 K (kO, X 1 - « y/iWx{x + 2)]- 1 / 2 J (A;O . 

Here Iq and K Q are the modified Bessel functions of the first and second kind, respectively. 
By continuing the perturbation analysis to higher orders we find that the approximate 
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solution x l - a t any (fixed) order has the following asymptotic form (as k — > oo) 

X z _ «F_(fc,Oe fef + F + (fc,Oe" fcf . (A5) 

For £ 7^ the functions F±(k, £) vanish at the limit /c — > oo. At the domain where both WKB 
and Green-Liouville approximations are valid, we equate Eq. (jA5j) with the corresponding 
WKB approximation [i.e., with Eq. (J23|) . extended to the required order]. Using this 
equation (for an arbitrary fixed order) we find that the term £a e - 25, o-25 2 approaches zero 
faster then any negative power of k, as k — > oo. We therefore conclude that the last term 
in the last brackets in Eq. (jAljl can be neglected in our calculation. Note, that the above 
mentioned equation can only provide us with a bound (which is sufficient for our purpose). 
This equation, however, can not be used to determine the value of the coefficient C2 (though 
it can be used to determine C\ up to a given order). The difficulty with the coefficient c<i is 
that it multiplies a subdominant term - a term which vanishes exponentially. 

We comment here that the entire perturbation analysis can be done with Green-Liouville 
approximation, this method has the advantage of being valid in the entire r > 2M region. 
However, we find that the WKB approximation is simpler, especially for the higher orders 
which are required by our analysis. 
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